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! Abstract 

In this paper we investigate the properties of astrometric microlensing of distant radio sources (QSOs 
' and radio galaxies) due to MACHOs, and discuss their implications for VERA (VLBI Exploration of Radio 
Astrometry). First we show that in case of astrometric microlensing of distant sources, the event duration is 
only a function of the lens mass and tangential velocity, but independent of the lens distance, in contrast to 
the well-known three-fold degeneracy for photometric microlensing. Moreover, the lens mass M is scaled by 
the tangential velocity v± as M cx v±, rather than M oc v\ which is the case for photometric microlensing. 
Thus, in astrometric microlensing the dependence of the lens mass on the unknown parameter v± is weaker, 
' indicating that the duration of astrometric microlensing event is a better quantity to study the mass of 
lensing objects than that of photometric microlensing. We also calculate the optical depth and event 
rate, and show that within 20° of the galactic center a typical event rate for 10 /ias-level shift is larger 
than 2.5 x 10~ 4 event per year, assuming that a quarter of the halo is made up with MACHOs. This 
indicates that if one monitors a few hundred sources for ^20 years, such an astrometric microlensing event 
is detectable. Since a typical event duration is found to be fairly long (5 to 15 years), the frequency of 
the monitoring observation can be relatively low, i.e., once per six months, which is rather reasonable for 
practical observations. We discuss practical strategy for observing astrometric microlensing with VERA, 
Qh' an d argue that an astrometric microlensing event due to MACHOs can be detected by VERA within a few 
q \ decades. 

Key words: Astrometry — The Galaxy: halo — Gravitational Microlensing — Radio Continuum: 
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that from astrometric observation one can determine the 
lens parameters such as lens distance and lens mass (e.g., 
Hog et al. 1995; Miyamoto & Yoshii 1995; Walker 1995). 
These astrometric observations are to be performed af- 
ter the event detection based on the photometric moni- 
toring, and in this sense, supplementary to photometric 
microlensing event. 

On the other hand, it is also known that astromet- 
ric observation alone can be used to study the nature of 
lcnsing objects. For instance, it has been argued that 
an astrometric effect of gravitational microlensing can 
be used to measure the mass of a nearby stars (e.g., 
Hosokawa et al. 1993; Paczynski 1996; 1998; Miralda- 
Escude 1996). Recently, the fluctuation of the radio 
reference frame by microlensing due to disk stars and 
MACHOs in the Galaxy has been studied by Hosokawa 
et al.(1997), and possibility to observe astrometric mi- 
crolensing events with ground-based or space optical- 
interferometer have been discussed by Miralda-Escude 
(1996), Bodcn et al.(1998), Paczynski (1998), and Do- 
minik & Sahu (2000). These studies revealed that as- 
trometric microlensing events are practically detectable 
and useful for studying the nature of lensing objects, if 
we obtain a position accuracy of 10 fias level. Thus, the 
astrometric microlensing will soon become one of major 
tools to study lensing objects. 

Following the recent studies on astrometric microlens- 
ing, in this paper we investigate astrometric microlensing 
of distant sources (QSOs and radio sources) due to MA- 
CHOs, focusing on the possible application to the VERA 
project (VLBI Exploration of Radio Astrometry). The 
VERA project (see for detail Sasao 1996; Kawaguchi et 
al. 2000; Honma et al.2000), which aims at radio astrom- 
etry at 10 lias level, has been approved in 2000, and is 
anticipated to start its observation by 2003. Distant ra- 
dio sources like QSOs and radio galaxies, which are one of 
major observational objects of VERA, are potential tar- 
gets for tracing astrometric microlensing with VERA, but 
few studies have been made concerning the astrometric 
microlensing of distant radio sources after the pioneer- 
ing work by Hosokawa et al.(1997). For these reasons, 
we study the astrometric microlensing of distant radio 
sources and discuss the implications for VERA. 

The plan of this paper is as follows; in section 2, we 
briefly review on astrometric microlensing, in particular 
parameters that can be extracted from the image trajec- 
tory. In section 3 we investigate the event duration of 
astrometric microlensing. In section 4 and 5, we calcu- 
late optical depth and event rate, and show that one can 
detect an astrometric microlensing based on monitoring 
of a few hundred sources for 10 to 20 years. Finally in 
section 6, we discuss the implication of astrometric mi- 
crolensing for VERA, and present a practical strategy for 
detecting astrometric microlensing. 



2. Brief Review on Astrometric Microlensing 

2.1. Astrometric Microlensing Event 

In this chapter we briefly review the characteristics 
of astrometric microlensing. When a source is closely 
aligned with a lens, the gravitational deflection of light 
produces two images whose positions are given by (e.g., 
Paczynski 1986) 



ri,2 = 



; ± Vu 2 + 4 



(1) 



Here ri,2 is the separation between the lens and the two 
images in the lens plane, and u is the separation between 
the source and the lens in the lens plane. Note that 
both ri t 2 and u are normalized with the Einstein ring 
radius Re; for instance, the lens-source separation u can 
be written as 

(2) 



u = 



Re ' 



where R is the lens-source distance projected onto the 
lens plane. As widely-known, the Einstein radius Re is 
given by (e.g., Paczynski 1986) 



Re — 



AGM D d D ds 



(3) 



Here M is the lens mass, and D d , D s , and D ds , are 
lens distance, source distance and lens-source distance, 
respectively. 

Here we introduce a coordinate system on the lens 
plane where the source is at the origin and the lens is 
moving parallel to the x axis (figure 1). The origin of 
time t = is set at the time of the closest approach 
between the lens and the source. In such a coordinate 
system, the lens position is described as (v±t, —/3Re), 
where v± is the tangential velocity of the lens, and (3 is 
the impact parameter normalized with i?E (note that (3 is 
always positive). With the parameters described above, 
the lens-source distance R can be written as 



R = V(v±t) 2 + (W 



(4) 



Here we calculate the positional shift of sources due to 
gravitational microlensing. Generally, two images pro- 
duced by gravitational microlensing cannot be resolved, 
but the centroid of two images can be observed to trace 
the astrometric effect in microlensing. If the lens is far 
from the source in the lens plane (i.e., u 3> 1), the sec- 
ondary image becomes too faint, and the image centroid 
coincides with that of the primary image. Since most of 
astrometric microlensing events with 10 /xas-level posi- 
tional shift have u 1 (as we will see in next sections), 
in the present paper we concentrate on the case of u 3> 1, 
and approximate the image centroid motion by that of 
the primary image (for more general treatment, see Bo- 
den et al.1998; Dominik & Sahu 2000). 
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Fig. 1.. Schematic view of the location of lens, source, and lensed 
image in the lens plane. Trajectory of an image is a circle in 
case of « > 1. Dots along the trajectory circle show the image 
position with a constant time interval, 0.4t q . 



In case that u 3> 1, the position of primary image 
(eq.p]]) is well approximated by, 



1 

u 



(5) 



The first term u corresponds to the lens-source separa- 
tion in the lens plane, and the second term l/u describes 
the image shift due to gravitational microlensing. Hence, 
the angular positional shift of the image relative to the 
original source position is obtained as 



u 



(6) 



Here 9e is the angular Einstein radius which is related to 
the Einstein ring radius Re as 



Re 
Dd' 



(7) 



The angular position of the image in the lens plane can 
be described as 



(0x,0 y ) = {O s cosS, O s smS), 
where the angle 6 is defined as (see figure 1) 
[3Re 



tan 8 



v±t 



(8) 



(9) 



(11) 



Using equations described above, one may rewrite the 
image position as 

-(6 E /m/ptE) qe/p \ 
(t/pt E y + i ' (v/w + iy ' 1 ' 

where we have introduced the Einstein ring crossing time 
4e which is defined as 

, _ Re 
t E = • 

Vj_ 

After a simple calculation, one may easily find the fol- 
lowing relation between 9 X and 9 Vl 

§H§)\ 

Hence, the trajectory of the image is a circle centered on 
(0, 6> E /2/3) with a diameter of 6 E //3 (Boden et al.1998; 
Dominik & Sahu 2000). We hereafter refer to this image 
trajectory as 'trajectory circle'. In figure 1, we show a 
schematic view of the image motion along the trajectory 
circle. The time intervals between adjacent dots are equal 
to 0.4x(/3£e). As seen in figure 1, the image motion 
is faster in the upper half of the trajectory circle than 
in the lower half, and vanishes toward the origin where 
the source should be located in case of no microlensing. 
The largest image shift parallel to the x-axis occurs when 
t q = ±/3£e; i- e -> R — n/2/3-Re- Note that an image goes 
through an upper quarter of the trajectory circle within 
the time interval of f q . Thus, i q gives a characteristic 
time-scale of an astrometric microlensing event. Note 
that from the observation of trajectory circle, one can 
obtain t q (— (3t E ) and 9e/ P- 

2.2. Angular Size of Astrometric Lens 

Here we estimate the angular size of astrometric lenses. 
Since we mainly consider the implications of astrometric 
microlensing for VERA, in the following analysis we as- 
sume that the sources are radio galaxies and QSOs which 
are much more distant (at least a few Mpc) than the lens 
(typically ~10 kpc). In such a case, the source distance 
D s and the lens-source distance D^s are much greater 
than the lens distance D^, and the Einstein ring radius 
(equation |3)) is well approximated as 



Re 



4GM 



D r] 



(13) 



Similarly, the angular Einstein ring radius 8 E may be 
approximated as 



AGM 



(14) 



If one substitutes typical parameters of MACHOs (e.g., 
Alcock et al. 1997; 2000), the angular Einstein ring can 
be estimated as 

^ n - i; — ) Ukpy mas - (is) 



0.5M, 



4 



[Vol. 53, 




Fig. 2.. Schematic view of the motion of source relative to lens. 
The large circle gives the lens size, which is equal to /9 max RE- 



Since the angular diameter of lens trajectory is given by 
0e/P, the maximum value of (3 for a detectable event is 
given by 



e/ #min; 



(16) 



where 9 m i n is the minimum position shift that can be 
measured. The astrometric projects planned to start 
within a decade (e.g., SIM, GAIA, VERA) aims at the 
position accuracy of 10 /Ltas or higher. Since we consider 
the observation of astrometric microlensing with those 
projects (in particular VERA), we set 8 m i n = 10 /ias. 
Thus, the maximum value of (3 (/3 m ax) is given as fol- 
lows, 



0n 



= 72 x 



M 



O.bMr. 



1/2 



-1/2 



8kpc J \ lOfxas 



(17) 



This indicates that the astrometric microlensing is de- 
tectable even when the minimum lens-source separation 
is larger than the Einstein ring radius by several tens. 
This makes the probability of astrometric microlensing 
much higher than the normal microlensing. 

3. Event Duration 

3.1. Basic Equations 

In this section, we investigate the event duration of 
astrometric microlensing. Usually, the duration of mi- 
crolensing event is defined as the time during which the 



source is inside the lens. For the photometric microlens- 
ing, the lens size is equal to the Einstein ring radius Re, 
and the event duration t p h is given as the time for which 
R < Re- Similarly to the photometric microlensing case, 
the lens size of an astrometric microlensing event is given 
by Pui^xRe, and the duration of an astrometric microlens- 
ing event is defined as the time for which R < /3 max RE 
(in other words, the time during which the image position 
shift is larger than the accuracy of position measurement 
0mm)- Thus, the event duration £ ast can be written as 

, 2/3 max i?E , . 

iast = COS0. (18) 

See figure 2 for a schematic view of lens-source geometry 
and the definition of the angle <f>. The average event du- 
ration can be obtained by integrating i ast over possible 
range of <$>. Since the probability of an event with a cer- 
tain value of cj> is proportional to cos <fi, the average event 
duration can be given as follows, 

/ iast cos <f> d<j> 



^a.st — 



J COS( 



Substituting equation ( |l8| ) into equation 
forming the integration with respect to q 
that 

7T (P n 



<Rf, 



v± 



(19) 

]l9| ) and per- 
, one obtains 



(20) 
(equations 



Substituting the expressions of Re and f3 B 
|p^[| , [ fl4"|| , and (l6)), one can obtain 

_ _ 2-kGM 

Interestingly, the average duration f ast contains only two 
unknown parameters, the lens mass M and its tangential 
velocity v±, but independent of the lens distance 
Hence, the degree of degeneracy in the event duration 
is lower than that in case of photometric microlensing, 
where the event duration depends on three parameter; 
v±, Z?d and M (note that this independence of the lens 
distance comes from the assumption of distant source. 
If this assumption is invalid, the event duration f ast de- 
pends also on the lens distance D^). Moreover, the av- 
erage event duration t ast is proportional to the lens mass 
M, and inversely proportional to the tangential veloc- 
ity v±. Thus, once the average event duration is ob- 
served, the mass is scaled as M cx wj_f as t. This de- 
pendence contrasts to that of photometric microlensing, 
where the event duration f p h is proportional to M 1 / 2 vJ 1 
(e.g., Paczynski 1986). Hence, for photometric microlens- 
ing events, the lens mass M is scaled asMa (wj_i p h) 2 , 
depending more strongly on the unknown parameter v± 
than in case of astrometric microlensing. This indicates 
that the event duration of astrometric microlensing is a 
better quantity to investigate the mass and the nature of 
MACHOs rather than that of photometric microlensing. 
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3.2. Halo Microlensing 

Here we estimate a typical value of the event duration 
for halo MACHO case based on the same parameters 
used in the previous sections. Substituting M — O.5M , 
v±_ = 200 km s -1 , and 9 m i n = 10 /ias, the event duration 
is obtained as 



15.3 x 



M 



O.5M 



/ V± \ 1 / flmin \ 

V 200km s"V V 10 Mas/ 



yr. (22) 



As seen in above equation, the event duration of astro- 
metric microlensing event is significantly longer than that 
of photometric microlensing event. This is due to large 
Anax for astrometric microlensing, typically ^70, while 
Anax is unity for photometric microlensing. In fact, if 
one takes a typical value of = 40 days from LMC 
microlensing (e.g. Alcock et al. 1997; 2000), a rough 
estimate gives that i as t ~ 2f3 max tE = 15 yr. 

An event duration of ~ 15 yr may seem to be too long 
for practical monitoring observation. This is because the 
event duration t ast is the time during which the image 
shift is larger than 9 m i n . On the other hand, we have 
already shown in figure 1 that the image trajectory is a 
circle, and that the image motion is faster in the upper 
half of the trajectory circle than in the lower half. As 
seen in previous section, tracing a quarter of the trajec- 
tory circle is enough to extract physical parameters for 
astrometric microlensing such as i q (= A^e) and #e/A 
Thus, practically it is necessary to trace the image mo- 
tion for the time interval t q , which was defined in section 
2 as the time during which the image passes the upper 
quarter of the trajectory circle. This duration t q is given 
by (see figure 2) 



t q = 0t E = 



Anax -R e i 



sm c 



(23) 



Similarly to t ast (eq.|nj), by performing the integration 
with respect to <j> we can obtain the average of t q as 



to = 



1/3 max Re 



(24) 



Comparing equations ([20|) and (|24|), one finds that t q — 
fast/"" j an d thus a typical value of i q is ~ 5 yr. This means 
that one can obtain physical parameters of astrometric 
microlensing such as A^e (= t q ) and 9e/(3 based on ~ 5 
yr monitoring of the event, which is rather reasonable for 
practical observations. 



3.3. Disk Microlensing 

Here we estimate the event duration for disk star lens- 
ing. As a typical mass of disk stars, we take the mass 
of M = 0.3M Q , corresponding to a lower main-sequence 



star. Since we mainly consider disk-star lensing at rel- 
atively high latitude (i.e., b > 10°) in which the lens is 
likely to be a nearby star, we take the tangential velocity 
of 20 km s _1 , corresponding to a random velocity of disk 
star in the vicinity of the Sun. With these parameters, 
the event duration for disk star microlensing is obtained 
as 



92 



M 



0.3M 



20km s" 1 / \10pas) 



yr. (25) 



Thus, the event duration for disk stars is much longer 
than that for halo MACHOs. This applies for most of 
the sky regions except for the Galactic plane, and hence 
the discrimination of disk lensing from MACHO lensing 
is easy. However, if one observes sources in the Galactic 
plane, the tangential velocity of the disk star is domi- 
nated by the galactic rotation, which is of ~ 200 km s _1 . 
This indicates that in the galactic plane the event dura- 
tion of disk lensing can be quite close to that of MACHO 
lensing, and hence sources in the middle of the galactic 
plane are not suitable for studying astrometric microlens- 
ing due to MACHOs, 

4. Optical Depth 

4-.1. Basic Equations 

Similarly to the optical depth of photometric mi- 
crolensing (e.g., Paczynski 1986), the optical depth for 
the astrometric microlensing can be defined as follows 
(Miralda-Escude 1996), 



ttA 



p(D d ) 



max^E ~n ^ d 



(26) 



Here p is the density, is the lens distance from the 
observer, and M is the lens mass. For simplicity here 
we assume that the lens mass M is unique, as usually 
assumed in microlensing studies. 

Substituting equations jl^), (|l4| ) and (|l6|), the optical 
depth (for distant sources) can be written as, 



7"ast 



16vrG 2 Af 



P dD d 



(27) 



As is already known (e.g., Miralda-Escude 1996), the op- 
tical depth for astrometric microlensing is proportional 
to the lens mass M, in contrast to that for photomet- 
ric lensing which is independent of the lens mass M. It 
is also remarkable that the optical depth is simply pro- 
portional to the column density, J p dD^. This simple 
dependence on the column density comes from the ap- 
proximation of distant source in equations (|l3| ) and (|l4|). 
In more general case, the optical depth is proportional 
to J p(l — x 2 ) dDd, where x = D d /D s is the fractional 
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Fig. 3.. Optical depth for astrometric microlensing due to MA- 
CHOs. The mass of MACHO is assumed to be O.5M . Thick 
line corresponds to (6, q)=(0°, 1), thin line to (0°, 0.75), thick 
dashed line to (10°, 1), and thin dashed line (10°, 0.75). The 
optical depth becomes larger with more flattened halo, because 
the halo density at the Galaxy's center increases with decreas- 
ing q. 




distance from the Sun (kpc) 

Fig. 4.. Differential optical depth for astrometric microlensing due 
to MACHOs. Notations are the same to those in figure 3. 



distance to the lens normalized with source distance (Do- 
minik & Sahu 2000). 

4-2. Halo Microlensing 

In order to calculate the optical depth due to MACHOs 
in the Galactic halo, here we introduce an axisymmct- 
ric logarithmic halo model (Binney & Tremaine 1987) to 
describe the halo density distribution. The density dis- 
tribution of the axisymmetric logarithmic halo model is 
given by, 

p(R,z) 



(2q 2 + 1)R 2 C +R 2 + (2-q- 



(R 2 + R 2 



z 2 q~ 
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Here R and z describes the cylindrical coordinate, and 
Vq, R c , and q are halo parameters which describe the 
circular rotation velocity, the core radius, and the halo 
flatness, respectively. In the following analysis, we set 
i? c = 5 kpc and vq — 200 km s _1 , which are typical 
values used in microlensing studies. We also introduce a 
cut-off radius beyond which the density vanishes. In the 
present study the cut-off radius is set to be 100 kpc, but 
note that our result would be only slightly changed by 
varying the cut-off radius, as we will see later. In addition 
to halo parameters in equation (^8|) , the mass of lens is 
necessary for optical depth calculation, as r ast depends 
on the lens mass (equation (27J). We set the MACHO 
mass M = 0.5M Q , according to the recent result from 
the LMC microlensing (e.g., Alcock et al. 1997; 2000). 
We also assume that the distance to the Galactic Center, 
i?o, is 8.5 kpc. 

The optical depth for astrometric microlensing is cal- 
culated for two values of q; q = 1 (spherically symmetric 
halo) and 0.75 (flattened halo). Figure 3 shows the opti- 
cal depth variation with Galactic longitude I. As seen in 
figure 3, the optical depth T ast is quite high toward the 
galactic center, being 2 ~ 3 x 10~ 2 . For a comparison, wc 
note that the optical depth for photometric microlensing 
T ph is 5.4 x 10~ 6 for / = 0° and b = 0° (q = 1 is assumed). 
This high optical depth for astrometric microlensing can 
be easily understood from the fact that typical /3 max is 
more than 70 (see previous section) . From the high value 
of T ast , one can expect that one of 30 to 50 distant sources 
seen around the galactic center is always being lensed, if 
the halo fully consists of MACHOs. Recent results of 
LMC microlensing indicate that the halo MACHO frac- 
tion is around 25~50%, and thus the optical depth may 
be 2 to 4 times smaller than what is obtained above. 
However, the probability is still as high as 5 x 10~ 3 , indi- 
cating that at least one of a few hundred sources is always 
being lensed. Note that in figure 3 the optical depth is 
larger for flat halo case (q = 0.75) than for spherical case 
(q = 1). This is because the halo density at the Galaxy's 
center increases with flattening of the halo. For instance, 
at the Galaxy center (R = 0, z = 0) the halo density is 
proportional to (2q 2 + l)/q 2 , which is 3.778 for q = 0.75 
while 3 for q = 1. 

In order to study a typical lens distance, in figure 4 
we plot the differential optical depth for astrometric mi- 
crolensing which is given by 



16ttG 2 M 
c ig 2 . ' 

u "min 



(29) 



inGq 2 



Figure 4 shows that the largest contribution to the optical 
depth occurs at D w 8.5 kpc, i.e., around the Galactic 
center. This is of course because that the halo density is 
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Fig. 5.. Optical depth for astrometric microlensing due to disk 
stars. The mass of disk stars is assumed to be 0.3Mq. From 
the top to the bottom, four lines correspond to b = 0° , 5° , 10° , 
and 20°. 
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distance from the Sun (kpc) 

Fig. 6.. Differential optical depth for astrometric microlensing due 
disk stars. Notations are the same to those in figure 5. 

largest at the Galactic center region. As seen in figure 4, 
the differential optical depth beyond 50 kpc is quite small 
compared to the Galactic center region. Therefore, the 
optical depth r as t only slightly depends on the value of 
cut-off radius (which is assumed to be 100 kpc), provided 
it is sufficiently large. 



4-3. Disk Microlensing 

Since the stellar disk is a dominant component in the 
inner Galaxy, we have to take into account the contribu- 
tion of stellar disk to astrometric microlensing. In order 
to calculate the optical depth due to disk stars, here we 
assume a disk which has an exponential profile both in 
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radial and vertical directions, namely, 

, p v ( R-Ro \z\\ 

Pd{R, z) — po exp I —I. (30) 

Here po is the local disk density, and d and h are radial 
and vertical scale length, respectively. We take typical 
values for these parameters as po — 0.08M Q , d = 3.5 
kpc, and h = 300 pc, which are the same to those used 
in previous studies (e.g., Hosokawa et al. 1997; Dominik 
& Sahu 2000). We show in figure 5 the optical depth 
for astrometric microlensing due to disk stars assuming 
M = 0.3M©. As seen in figure 5, the optical depth is 
close to 0.1 at the galactic center, indicating that one 
of about 10 sources is always being lensed. This optical 
depth by disk stars is consistent with previous studies 
such as Hosokawa et al.(1997). If one observes a source 
relatively far from the galactic center but in the middle of 
the galactic plane (6 = 0°), the optical depth due to disk 
star is larger than 0.02 for I < 50° . Therefore, as long 
as one observes sources in the Galactic plane, the optical 
depth is significantly higher than that for halo MACHOs, 
and hence sources in the galactic plane are not suitable 
for astrometric microlensing study of MACHOs (but of 
course they are useful for studying the galactic disk, par- 
ticularly the stellar mass and the disk scale length; e.g., 
Hosokawa et al. 1993; Dominik & Sahu 2000). How- 
ever, as seen in figure 5, the optical depth due to disk 
stars decreases quite rapidly with increasing the galactic 
longitude. Thus a source at relatively high galactic lat- 
itude (b > 10°) may be used to search for astrometric 
microlensing due to MACHOs. 

Similarly to the halo case, in figure 6 we show the dif- 
ferential optical depth for disk stars. As easily under- 
stood, the differential optical depth in the galactic plane 
(b = 0°) becomes maximum at the Galactic center. On 
the other hand, for sources at relatively high galactic 
latitude, the differential optical depth for disk stars is 
largest at D<j = 0, and hence it is dominated by local 
stars located within a few kpc from the sun. In such a 
case, one can discriminate the disk star microlensing from 
halo MACHO microlensing based on the event duration, 
as seen in the previous section. 

5. Event Rate 

The event rate is another important quantity in mi- 
crolensing study, because the event rate is tightly re- 
lated to the observational strategy of microlensing search. 
Once the optical depth r ast and the average event dura- 
tion fast are known, the event rate r ast is obtained as, 

r ast = (si) 

"ast 

Substituting equations (|Tj) and ( p7| ) into above equa- 
tion, the event rate for astrometric microlensing may be 
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written as follows, 

r as t = ^^ f P dD d . (32) 

Remarkably, the event rate r as t does not depend on the 
MACHO mass M, as is the case for the optical depth for 
photometric microlensing (r p h). Also, the event rate r as t 
is simply proportional to the column density J p dD^. 
These facts indicate that the event rate of astrometric mi- 
crolensing is a good tracer of density distribution rather 
than the mass of individual lenses. 

If we put into equation (|3l| ) F as t = 15 yr and T ast = 
1.5 x 10~ 2 , which corresponds to the optical depth at 
20° away from the Galactic center, we obtain 

T ast = 1.0 x 1(T 3 f — T ) 

\200kms _ V 

x (^-) 1 event yr- 1 . (33) 
\10 pas ) 

Thus, if one observes 100 sources for 10 years, at least 
one event will be detected. However, note that this cal- 
culation is for the case that the halo totally consists of 
MACHOs. Since the halo MACHO fraction is estimated 
to be 25%-50% by recent studies (Alcock et al. 1997; 
2000), one has to observe 2-4 times more sources or longer 
period to detect an event. This indicates that it is neces- 
sary to monitor a few hundred sources for a few decades. 
This number of source and the time span may seem to 
be long, but the frequency of monitoring can be fairly 
low because of the long event duration, and thus such 
an observation can be done with upcoming astrometric 
projects like VERA (we will discuss this in next section 
in detail). 

For the disk star lensing, the event rate is expected to 
be much smaller than that for halo MACHO lensing be- 
cause the event duration i ast is significantly larger. From 
figure 5, the optical depth for disk star lensing at b = 20° 
is 1.2 x 10~ 3 (note that this value is almost independent 
of the galactic latitude I). Using the event duration of 92 
yr, one obtains an event rate for astrometric microlensing 
due to disk star as 1.3 x 10 -5 , which is smaller by two 
orders of magnitude than that for halo MACHO lensing. 
Thus, the contribution of disk star lensing is negligible in 
the off-plane region. 

6. Implications for VERA 

In this section we discuss the implication of astromet- 
ric microlensing for VERA. As we have briefly mentioned 
in section 1 , VERA is a new VLBI array to measure the 
relative position of two adjacent sources based on the 
phase-referencing technique. Although VERA aims at 
astrometry with 10 pas- level accuracy similarly to other 
space missions (e.g., SIM, GAIA), VERA is very unique 



as it is an array of ground-based radio telescopes. The ad- 
vantage of using VERA for astrometric search is that the 
project life time is much longer than that of other space 
missions, being ~20 years (as for space missions like SIM 
and GAIA, the mission lifetime is around 5 years). As 
we have seen in the previous sections, the event dura- 
tion is as long as 10 years, and thus a monitoring span 
exceeding 10 years is indispensable for studying astro- 
metric microlensing event. 

As we have shown previously, to detect an astromet- 
ric microlensing event, a few hundred sources should be 
monitored for a few decades. If one wants to detect more 
events for statistical studies, the number of sources to 
be observed becomes larger. Thus, it is important to 
examine if there are enough number of sources for astro- 
metric microlensing search with VERA. The main tar- 
get of VERA are compact continuum sources and maser 
sources. The compact continuum sources are mainly dis- 
tant radio galaxies and QSOs, and thus they are suit- 
able to astrometric microlensing studies. On the other 
hand, the maser sources are mostly late-type stars and 
star forming regions in the Galaxy, and so unsuitable 
for the purpose of the present study. As for continuum 
sources, about 2000 VLBI sources (which are compact 
enough to observe with VLBI) are currently known (e.g., 
ICRF catalog, Ma et al. 1998; VLBA Calibrator Survey, 
Peck & Beasley 1998). At high galactic latitude where 
VLBI sources are searched relatively deeply, the average 
number of VLBI sources is close to 0.08 per square degree. 
This number roughly corresponds to one VLBI source in 
a circle with radius of 2°. If the same source distribu- 
tion applies to the sky region toward the galactic center, 
one can expect to find about 100 sources within 20° of 
the Galactic center, and 220 sources within 30° of the 
Galactic center, indicating that there are enough num- 
ber of VLBI sources that can be observed with VERA. 
Moreover, there are several hundreds of VLBI source can- 
didates like those in Jodrel-Bank VLA Astrometric Sur- 
vey (JVAS; Patnaik et al. 1992; Browne et al. 1998; 
Wilkinson et al. 1998). By a deeper and more massive 
survey, one can expect to increase the number of VLBI 
sources that can be used for astrometric microlensing 
search. Hence, one can expect to have a plenty of VLBI 
sources around the galactic center regions. 

Next we estimate the time required for such a moni- 
toring observation. The signal-to- noise ration (S/N) of 
VLBI can be estimated using the following equation (e.g., 
Thompson et al 1986), 

fls/N = ^7\ L Al i 2 U S. (34) 

Here A is the effective antenna collecting area, T sys is 
the system noise temperature, Ai/ is the bandwidth, t\ is 
the integration time, k is the Boltzmann constant, r\ is 
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the quantization efficiency, and S is the source flux (sub- 
scripts denote the observational stations). The major 
specifications of VERA telescopes are currently antici- 
pated as follows; the effective collecting area A of 157 m 2 
(assuming 20m diameter and antenna efficiency of 0.5), 
system noise temperature T sys of 120 K, bandwidth Av of 
256 MHz, and the quantization efficiency r\ of 0.88 (two- 
bit quantization). With these values, the S/N ratio Rs/n 
is <~ 16 for 100 mJy source with integration time U of 300 
seconds. Since the accuracy of fringe phase measurement 
is determined as l/i?g/N m radian, its accuracy will be 
0.063 radian, or 3.6 degree. This corresponds to an error 
in geometric delay of 0.13 mm (— 3.6 x 7r/180 x A, where 
A is the wavelength) at the wavelength of 1.3 cm, which is 
one of the major frequency band of VERA. The position 
accuracy can be roughly estimated as (0.13 mm/2300 
km)=5.65x 10~ n radian, which equals to 11.7 /xas (2300 
km is the longest baseline of VERA). Thus, for sources 
with a flux of ~ 100 mJy, the thermal error can be re- 
duced to 10 fxas level within 5 minutes. 

If one spends 20 minutes for a single measurement (one 
measurement of separation of two adjacent sources), the 
total number of sources that can be observed in 24 hours 
exceeds 100 (note that two sources are observed in a sin- 
gle measurement). Thus, 2-3 days are enough for observ- 
ing a few hundred sources. In practice, one separation 
measurement for each pair is not enough for detecting 
a position shift of ~10 ^as. To achieve higher position 
accuracy (i.e., a few /ias), we may need 10 sets of ob- 
servations of a few hundred sources, which requires total 
observation of 20-30 days. On the other hand, typical 
timcscale for a image to move along a quarter of the 
trajectory is about 5 years, and so the interval of such 
observations can be as long as six months. If we spend 
30 days for astrometric microlcnsing in six months, this 
requires 15% of the telescope time to be dedicated to 
such a monitoring observation. This is indeed massive 
expense, but is worth spending when one considers the 
significance of scientific outputs; in addition to astromet- 
ric microlensing, one can establish a radio reference frame 
with an accuracy of 10 fxas, which is more accurate by two 
orders of magnitude than the reference frame currently 
achieved, and also one can measure proper motions of 
nearby AGNs (up to a few tens Mpc), from which 3 di- 
mensional structure of the local universe can be studied. 

One possible way to incorporate such a massive mi- 
crolensing search with VERA is to perform a geodetic 
observation based on the dual-beam observation. The 
geodetic observation is necessary to determine the base- 
line parameters with sufficient accuracy. Usually, a 
geodetic observation is made at S (2GHz) and X (8GHz) 
bands with a single beam. In case of S/X geodetic ob- 
servation, the baseline parameters and absolute position 
of radio sources with an accuracy of 1 mas are solved 
based on a few hundred observations of radio sources. 



On the other hand, baseline measurements are also pos- 
sible based on the relative fringe phase of two sources that 
are observed with dual beam system of VERA. In case 
of dual-beam geodetic observation, what can be obtained 
are the baseline parameters and 'relative' position of ad- 
jacent two sources with an accuracy of 10 ^as level. Thus, 
if dual-beam geodetic observation is performed, one can 
obtain relative position of radio sources for a few hundred 
sources with an accuracy of 10 fias. Currently, VERA is 
planned to perform one-day geodetic observation every 
week, and thus the frequency of the observation is just 
the same to what we need for astrometric microlensing 
search (i.e., 15% of the telescope time). Therefore, dual- 
beam geodetic observation, if possible, is more practical 
for the massive monitoring of astrometric microlensing. 

7. Conclusion 

We have investigated the properties of astrometric mi- 
crolensing of distant sources due to MACHOs. We have 
shown that in case of astrometric microlensing of distant 
sources, the event duration depends only on two unknown 
parameters, the lens mass M and the tangential veloc- 
ity vj_, in contrast to that for photometric microlens- 
ing of nearby stars. We have also shown that a typical 
event duration of astrometric microlensing is 15 years, 
and also that the event rate of astrometric microlcnsing 
is larger than 2.5x 10~ 4 within 20° of the Galactic center, 
assuming that 25% of Galaxy's halo is made up with MA- 
CHOs. This implies that one can detect an astrometric 
microlensing event if a few hundred sources are monitored 
for 20 years. Detections of astrometric microlensing will 
be possible with VERA within decades, and thus astro- 
metric microlensing will become a new tool to study the 
nature of MACHOs. 
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